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Abstract 

We calculate the rate of decrease of the expectation value of the transverse com- 
ponent of spin for spin-1/2 particles in a magnetic field with a spatial gradient, to 
determine the conditions under which a previous classical description is valid. A den- 
sity matrix treatment is required for two reasons. The first arises because the particles 
initially are not in a pure state due to thermal motion. The second reason is that each 
particle interacts with the magnetic field and the other particles, with the latter taken 
to be via a 2-body central force. The equations for the 1-body Wigner distribution 
functions are written in a general manner, and the places where quantum mechanical 
effects can play a role are identified. One that may not have been considered previ- 
ously concerns the momentum associated with the magnetic field gradient, which is 
proportional to the time integral of the gradient. Its relative magnitude compared 
with the important momenta in the problem is a significant parameter, and if their 
ratio is not small some non-classical effects contribute to the solution. 

Assuming the field gradient is sufficiently small, and a number of other inequal- 
ities are satisfied involving the mean wavelength, range of the force, and the mean 
separation between particles, we solve the integro- partial differential equations for 
the Wigner functions to second order in the strength of the gradient. When the same 
reasoning is applied to a different problem with no field gradient, but having instead 
a gradient to the z-component of polarization, the connection with the diffusion co- 
efficient is established, and we find agreement with the classical result for the rate of 
decrease of the transverse component of magnetization. The corresponding result for 
this rate in the absence of collisions is much greater. 

An approximate value for the 3 He self diffusion coefficient based on a hard core 
potential is found to be 23.1 mm 2 /sec at room temperature and 7 atmospheres. This 
falls between the results obtained in an NMR experiment having these values of tem- 
perature and pressure, in which the spatial gradient of the field is time dependent, 
and earlier experiments using other techniques, corrected for temperature and pressure 
variation. 



1 Introduction 

The problem of a system possessing magnetization in the presence of a magnetic field 
having a spatial gradient to its strength was considered classically by Carr and Purcell 
[jij, Torrey P], and Abragam |J. The rate of precession of the transverse component of a 
particle's magnetic moment depends on where it is located in the field, but thermal motion 
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and collisions smear this out. This was taken into account in Q by carrying out a random 
walk, and in Q and Q by adding a diffusion term to the equation for the time rate of 
change of the magnetization. It leads to a reduction of its transverse component which 
is exponential in the cube of the time. Stejskal and Tanner B took the same approach 
but generalised the problem by permitting the gradient to be time dependent. In the 
approximate solution obtained there j| the reduction of the magnetization has a more 
complicated time dependence. A classical treatment of NMR experiments involving spatial 
variation of the magnetization, including the effect of RF pulses, is given in Sodickson and 
Cory |J. We shall examine the diminution of transverse magnetization as a quantum 
mechanics problem to determine under what conditions the classical approach described 
above is justified. 

An isolated spin 1/2 particle is described by a wave function the spin part of which 
can be written as a linear combination of two basis states, which can be taken to be 
eigenstates of the z-component of spin with eigenvalues ±1/2. Any given values of the 
relative magnitudes and phases of the coefficients of those two states represents a particle 
whose spin is pointing in some definite direction in 3-dimensional space. In the presence 
of a constant magnetic field the spin rotates coherently about the field. 

The theoretical description of the process of interest is more complex and proceeds in 
two steps, precisely following the order in which the experimental conditions are created. 
Focussing on a single particle, it is located in a heat bath, which means there is a classical 
probability distribution for its momentum ||. In addition, the method of preparation 
(partially polarized) creates a probability distribution for its spin direction 0]. At this 
first stage already the momentum and spin variables for the particle must be described 
by a density matrix @, 0, ^ • In the presence of a constant magnetic field the probability 
distribution of the particle's spin rotates coherently about the field direction, just as for 
the case of an isolated particle in a pure state. At this stage the central interaction of the 
particle with the other particles, which are in the same heat bath, has no effect on the 
spin. 

In the second step a gradient to the strength of the field is turned on. The momentum 
distribution of the particles becomes distorted away from thermal equilibrium because a 
magnetic moment in the direction of the field experiences a force in the direction of the 
gradient. Now the 2-body interaction does have an effect. The expectation value of the 
transverse component of the particle's spin is determined by the off-diagonal elements of 
the density matrix in spin space, and we shall see that the rate at which those elements 
decrease with time is affected by the interplay between the two-body interaction and the 
field gradient. 

The effect of (anti)symmetrization on transport properties was studied by Lhuillier and 
Laloe and included initially in the present work. Those papers were not concerned 
with the effects of an applied magnetic field; and there appears to be a technical difference 
between their treatment of the effect of 2-body scattering on the 1-body density matrix 
and that given in Snider jn]], which we follow. This is discussed in the text. 

In the Theory section we formulate the problem in terms of the 1-body Wigner distri- 
bution functions, and the coupled integro-partial differential equations that they satisfy 
are obtained as well as the expressions for the expectation value of the particle's spin in 
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terms of those functions. Some important parameters are identified, including the mo- 
mentum associated with the magnetic field gradient. Provided the gradient is sufficiently 
small and a number of other inequalities are satisfied, the equations are solved in Section 
3 by treating the gradient of the magnetic field as a perturbation, to second order in its 
strength. The successive terms of the solution are gathered together in Section 4, where 
the result is compared with the classical solution. The approximations encountered along 
the way are summarized in Section 5. 

A word about terminology. The bulk of the paper is concerned with the quantum 
mechanical calculation of the expectation value of the spin of a particle. The classical 
calculation with which it is compared speaks about magnetization, which is the magnetic 
moment per unit volume. The two quantities are simply related: multiplying the expecta- 
tion value of the spin by the magnetic moment of a particle and by the number of particles 
per unit volume gives the magnetization since all the particles are identical. 



2 Theory 



The Hamiltonian for the system is taken to be a sum of 1-body terms, plus a sum 

of 2-body potentials V. H^ 1 ' contains the kinetic energy plus the interaction energy of 
the particle's magnetic moment with the applied magnetic field, B z , taken to be in the 
z-direction. For spin-1/2 particles this becomes 



p 

2M 



7 



ha z B, 



(1) 



where M is the mass of the particles, 7 is the gyromagnetic ratio and a z is the Pauli spin 
matrix. With the gradient of the magnetic field strength in the direction of the unit vector 
u, B z can be written 



B z = B + u-xG(t). 



(2) 



Bq is the constant part of the field and G(t) represents the magnitude of the gradient of 
the field. 

We follow the development in Snider (llj , but with the inclusion of (anti)symmetrization. 
There the density matrix for an N particle system is reduced by taking successive traces 
down to the 1-body density matrix, p. It satisfies the equation 



dp(l) 
dt 



+ TV 

(2) 



V(l,2),pW (1,2) 



(3) 



where the arguments of the operators refer to particles number 1 and 2 and p^ is the 
2-body density matrix. (A superscript (1) is omitted on the 1-body p.) It is indicated 
in the second term within the parenthesis in (Q) that the trace is to be taken over the 
variables of particle 2. The hierarchy of equations is cut off by neglecting the terms in 
the equation for dp^ /dt which involve p( 3 ' , arguing that in a dilute gas when particles 1 
and 2 are interacting the interaction with a third particle can be neglected [10]. With n 



being the number of particles per unit volume and a the range of the two-body potential, 
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this requires that a <C d = n -1 / 3 . Normalising to the conditions found in the experiment 
reported by Schmidt et al. [^] involving a partially polarized gas of 3 He atoms, which 
was performed at room temperature and a pressure of 7 atmospheres, 

f^cisc-^^n ') ,4) 

d \7atm T J \2A-10- 10 mJ KJ 

Although the ratio a/d may be on the borderline of being sufficiently small for that partic- 
ular experiment |l2| , we shall neglect the possibility of triple collisions, and will continue 
to take a = 2.4- 10 -10 m as an estimate of the range of the He-He potential when examining 
the size of other relevant parameters. 

We assume that prior to a collision between particles 1 and 2 they are uncorrelated, 
and hence p( 2 ) is just the product of the two 1-body density matrices; this corresponds to 
the Boltzmann property for the classical distribution functions. Using scattering theory 
the effect of a collision can be obtained as 

pp>(i,2) = nWp(i)p(2)nMt (5) 

where is the wave operator [|ll], f[0| . When this expression for p^ is inserted into (|3[) 
it becomes an equation for the 1-body density matrix p alone. We shall seek a solution 
of this equation to second order in the strength of the gradient G(t). In Lhuillier and 
Laloe H instead of the wave matrix f^ + ) in Eq.@ the S- matrix is used. This relates the 
density matrix in the distant past, before the collision, to its value in the distant future 
after the collision is completed. But we are interested in the time development of p during 



a collision [1C] and choose the wave matrix, therefore. Our numerical results would be 
different had we chosen the S'-matrix. 

In the experiment reported by Schmidt et al. |jl2[ ; a gas of 3 He atoms is partially 
polarized in the z-direction and fairly uniform spatially. At t = it is tipped through an 
angle 9 thereby producing a component of spin in the plane transverse to the direction of 
the magnetic field as well as a component in the z-direction. It is the diminution of the 
transverse component with time that was studied, and which we wish to evaluate. The 
most general form of the density matrix for spin- 1/2 particles can be written 

p = pil + pcr 

= pil + p z a z + -(p + a- + p-a + ) (6) 

where for any operator, A± = A x ± iA y , and the are the Pauli spin matrices. The 
four operators pi, p z , p + , and p_ = p\_ operate only in ordinary (position or momentum) 
space. 

It is convenient to work with the Wigner distribution function || generalized to a 
2x2 matrix to include the spin degree of freedom. This generalization was also made in 
reference Q . It is defined by 

1 f d 3 q , 
/ mm /(x,p,t) = J < P ~ q> m \p\P + q,m / > exp(-2ix • q/fi) (7) 

= y rf 3 y < x - y,m|p|x + y,m' > exp(2ip • y//i) (8) 
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where the indices m and m' take the values ±1/2. It follows from the fact that p is 
hermitean that the 2x2 matrix / is also. We also need the inverse relation 

< p,m\p\p',m >= (2-Khf J d 3 xf mm ,(x, (p + p')/2, t) exp[i(p' - p) ■ x/ft]. (9) 

Taking the partial derivative of Eq.(|7|) or (H) with respect to time and using Eq.(|3]) 
gives 

mm' _ qK i 07 i qV /in', 

g t '-'mm' ~r 10 mm' ' mm' \ LyJ ) 

where S K arises from the kinetic energy term in , 5" 7 from the magnetic energy term 
in ifW ( a nd S v from the 2-body potential energy term. S K is the well known flow term 

Smm> = ~~ J^P " V/ mm /. (11) 

To obtain S* 7 requires evaluation of the commutator [a z B z , p\. Making use of Eqs. (|2|) and 
(|) yields 

S l/2 = -^G(t)u-V p / 1/2 (12) 

Sll/2 = ^G(t)u.V p /_ 1/2 (13) 

5j = -i 7 5z/+ (14) 

where we have introduced abbreviated notations for all the S terms as follows: SV/2 = 
5*1/2 1/2 ; -5-1/2 = 5-1/2 -1/2 ! an d 5'_|_ = 5_!/ 2 i/2! an d the same definitions apply to the 
components of the Wigner distribution functions /. We shall be particularly interested 
in the function / + , which involves the off-diagonal elements of p in the 2-dimensional 
spin space. In conjunction with Eq .(p!o[ ) for the time rate of change of the distribution 
functions, the interpretation of Eqs. (|12|) and ( |l3|) is as follows. A magnetic moment having 
a component in the direction of (or opposite to) a magnetic field whose strength has a 
spatial gradient will experience a force in the direction of (or opposite to) that gradient. 



Eq.(14) says that a magnetic moment transverse to the direction of the magnetic field 
undergoes Larmor precession about the field. 

With the density matrix normalized to Tr(p) = N, the total number of particles, it 
is straightforward to obtain the following results for n, the number of particles per unit 
volume, and < er >, the expectation value of the particle's spin at position x and time t 



n 



|d 3 p[A /2 (x,p)+/_ 1/2 (x,p)] (15) 



and 



f d 3 p[f 1/2 (x,p) - /_i/ 2 (x,p)] 
2 /rfM/i/2(x,p) + /-i/ 2 (x, P )] 1 j 

. 2/d 3 p/ + (x,p,t) 

+ /dWi/ 2 (x J p)+/-i/ 2 (x,p)r 1 1 
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Once these quantities have been obtained the net magnetic moment per unit volume, M, 
can be obtained simply by multiplying < a > by the magnetic moment of a single particle, 
jh/2, and the number of particles per unit volume, n. 

We now turn to S v , the collision term in Eq.([lO|). Introducing the spin and momentum 
variables for particle 2, over which the trace indicated in Eq.(||) will be taken, and the 
antisymmetrization operator (/ + eP\2) where P\2 permutes all the variables of particles 
1 and 2 and e = — 1 for fermions, S v can be written 

v , . i 1 x - f d 3 q d 3 p 2 



■ < p-q,m;p2,m 2 |(/ + ePi 2 )[^,p (2) (l, 2 )]|p + q ) ' m ' / ;P2,^2 > 

• exp(-2ix • q/h). (18) 

It is sufficient to antisymmetrize just one of the states because V and p^ are symmetric 
under the interchange of the two particles. In the Appendix we will develop the collision 
term making two approximations. The first is to assume that the 2-particle potential 
energy V consists purely of a central interaction. For example, at a separation of 2.4 A 
the moment-moment interaction energy of two 3 He atoms is ~ 10 -8 mev, whereas the 
central interaction is ~ 1 mev. The second approximation amounts to assuming that the 
gradient of the magnetic field strength is sufficiently small; this comes about as follows. 

Examination of Eqs.(|i~0|) through (|i~4|), which represent the contributions to the time 
derivative of the respective distribution functions, shows that only $7 contains dependence 
on position via the gradient term in the magnetic field, B z . Neglecting the collision terms 
S v for the moment, would lead to the following conclusions. With the particle density 
n and the expectation value of the spin initially independent of position, n and < a z > 
remain so since f±i/2 are independent of position, n is, of course, also independent of 
time, and < a z > as well since a z commutes with the Hamiltonian. On the other hand /+ 
would contain a factor exp[— i^{BQt + u • xi ? (i))] as the only position dependence in the 
problem, where 

F(t)= ( dt'G(t'). (19) 
Jo 

We show below that all these conclusions remain valid with collisions and write, therefore, 
/+(x,p,t) =0 + (p)exp[-t7(fl o * + u-xF(*))]. (20) 

The product jF(t) is the wave number associated with the magnetic field gradient ||; in 
that reference it is designated as k. 

The consequences for the matrix elements of p can be seen as follows. From Eq.(|9]) 
those matrix elements with m = m' , since they arise from f±i/2 which do not depend on 
position, are local in momentum space, 

< p,m\p\p',m >= (27rn) 3 / m (p)(2vr/i) 3 ^ 3 )(p' - p). (21) 
On the other hand, for the off-diagonal matrix elements we have 

< p, -l/2|p|p', 1/2 >= (2^) V((p + pO/2) exp(-^ot) 

•(2vr?i) 3 <5 (3) (p'-p-7?iuF(t)). (22) 
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Position dependence of a distribution function is sometimes neglected in the collision 
term on the assumption that the quantity that is varying (the strength of the magnetic 
field in the present case) changes little over the distance of a mean free path But 
the (5-function in Eq.(f22|) contains another relevant quantity, jnuF(t), that has nothing 
to do with the mean free path. It is a momentum that arises from the oppositely directed 
forces felt by particles with their moments in the direction of the field, or opposite to it, 
due to the field gradient. After inserting two complete sets of 2-body spin and momentum 
states into Eq.(|i~8|) matrix elements of the type < p, — l/2|p|p', 1/2 > appear linearly 
in the collision term and bilinearly (or not at all) in S±y 2 - When the 5- function 
from Eq.(^) is combined with the two momentum conserving 5-functions from the t- and 
fiW-matrix elements [see ( A-l ) and (|A-2| )1 it has the following effects. In Eq.(|i~8|) for 



<S*Xi/2 the integration variable q is required to be zero; but in it becomes jhF(t)u/2. 
Consequently acquires a factor exp[— i^(Bot + u • xF(t))], and as a result every term 
in Eq.(|l(]) for df+/dt has that factor. This confirms the statements above that f±i/2 are 
independent of position and /+ has the structure shown in Eq.(|20|). No assumption about 
the magnitude of F(t) has been made at this point. 

If the magnitude of jhF(t) were not sufficiently small compared to the important 
momenta in the problem there would be significant quantum mechanical consequences for 



the collision terms of the displaced argument of the 5- function in Eq.(22). The arguments 
of 5+ and the arguments of the i-matrix elements would be shifted by amounts of order 
7?ii ? (i)u. We shall delay discussion of the magnitude of the effect on the collision terms 
due to the altered t-matrix elements, but note here that it would even lead to the presence 
of non-energy-conserving matrix elements. 

The effect of the variation in g + can be seen directly. Although the distribution 
functions are of the Fermi-Dirac type, we shall be interested in conditions of negligible 
degeneracy, nX <C 1, where A is the mean wavelength of the particles. In that case we can 
use the Maxwell-Boltzmann distribution, and its momentum derivative is of order (p/p 2 ) 
times the function itself, where p is the mean momentum at temperature T. The fractional 
shift in the magnitude of g+, therefore, is of order "fhF(t)/p. We shall estimate it by 
normalising to the experimental conditions found in reference [jT^]. The gyromagnetic ratio 
for 3 He is 204 MHz/Tesla, and the maximum value of F{t) in that experiment, designated 
Fm, is 3.5 -10~ 7 Tesla-sec/cm. Estimating p by the mean value of the momentum at 
temperature T, we obtain 

25M K9 . 5 . 10 -»m(^)(^^) 1/! . (23) 
V \7>ho/ V F M A M T I ' 

For any other set of experimental parameters it is only necessary to insert the appropriate 
values of the gyromagnetic ratio 7, the temperature T, the molecular mass M, and the 
integrated strength of the gradient of the magnetic field, F(t). Note that a larger strength 



of the field gradient or a smaller temperature would increase the ratio in Eq.(23), but the 
variation with T is slow. 

For the remainder of the paper we shall assume that the shift in the arguments of 
g + and the t-matrix elements can be neglected, and this leads to the expressions for the 
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collision terms presented in Eqs. QA-17 ) and ( A-18| ). Eqs.(|K]) then become the following 
set of coupled nonlinear integro- partial differential equations. 



9/i/2 (p) 



Of. 



-^G(t)u.V p / 1/2 (p) + S| / /2 (p 
<9/-i/2(p) _ 7 



at 



ftG(i)u-V p /_ 1/2 (p) + Sr i/2 (p) 



(24) 



(25) 



and in the equation for df+/dt the derivative of the exponential with respect to time 
cancels the term Si leaving 



dg+(p) . 7 r?,,s , , . -c v r \ 



(26) 



3 Perturbation Expansion 

Our approach for solving these equations, consistent with the approximation described 
above, is to treat the strength of the field gradient, G(t), as a perturbation. Since it is 
multiplied by the gyromagnetic ratio 7 we shall expand the distribution functions up to 
second order in 7, 

/±i/2 = f% + 7/2/2 + l 2 ftl,2 + Oil') (27) 

and 

5 + = ^° ) +75l 1) +7 2 ^ 2) +0(7 3 )- (28) 

This procedure converts the collision terms S v from being bilinear in the functions / and 
g to being linear in the successive approximations /W and 



Zeroth Order 



The terms of order 7 in Eqs.(p4]) - (26) represent the initial condition in which the prob- 
ability of finding the particle with a given momentum is given by a thermal distribution, 
and the expectation value of its spin is specified; both quantities are initially independent 
of position. Those distribution functions can be written 

/|J /2 (p) = n ±1/2 B( M \p) (29) 

and 

g^(p) = n + B^\p), (30) 

where 

B M(p) = (/?( m Vtt) 3/2 eM~P m P 2 ) (31) 

is the the Maxwell-Boltzmann distribution for a particle of mass M, and/?( M ) = (l/2MkT). 
As expected from conservation of energy |l3}| , when these zeroth order distribution func- 
tions are inserted in the collision integrals they can be shown to vanish, confirming that 
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Eqs. (24) - (26) are satisfied to this order. The unitarity relation gets used here and in the 
subsequent development [see Appendix]. The significance of the normalization constants 
can be seen from Eqs.([T5D to (p"7|), since they determine the initial values of the number 
density, n, and the expectation value of the spin, < a >: n = n\i 2 + ra_i/2> < a z >= 
( n i/2 ~~ n -i/2)/ n ) an( i < a + >= n + /n. 

First Order 

At the next order in 7 the distribution functions are distorted in the direction u of the 
gradient, so we write 

<•(!) f n \ _„.... r(M)/ x . (i) 



/S /2 (p) = n ±1/2 B( M \p)u • ph% 2 (32) 



and 
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^(p) = n + B( M \p)u.ph^. (33) 

In principle the functions hP 1 ' could depend on the magnitude of the momentum, and a 
systematic expansion in a complete set of Sonine polynomials, for example |^], could be 
undertaken. We shall take the functions to be independent of p, thereby limiting that 
expansion to just the lowest term. Note that these first order terms in the distribution 
functions contribute nothing to the number density and the expectation value of the spin, 
Eqs. ( |l5|) - (17), because their integrals over all momenta vanish. As expected if < a z >^ 
the particles acquire a non-zero expectation value of momentum, and it can be shown to 
be proportional to the product of («i/2 — n_i/ 2 )/n an d jhuF(t). The latter factor is the 
same one we required to be small compared to thermal momenta in the discussion above. 



Now pick out of Eqs. ( |24| ) - (26) all the terms of first order in 7. Since the functions 
depend only on the time, these become coupled ordinary first order linear differential 
equations 

dh {l) 

n 1/2 B( M Hp)u ■ = n 1/2 np^B^(p)n • P G(t) + (34) 

dh (1) 

n_ 1/2 B( M \p)u ■ p—^ = -„^«#)(p)u • pG(t) + S V [% (35) 

and 

M (J °(p)n • = 4m M (p)u • pF(t) + S V + {1) . (36) 

The expressions for the first order collision terms S v ^ are given in the Appendix, where 
the above equations are developed into Eqs .( [A- 24 ) to ( |A-26 ). 



There are three quantities in those equations, and that depend upon the 

details of the 2-body potential energy; see Eq.( A-23] ). Iu, which is the average over a 



Maxwell-Boltzmann distribution of the cube of the momentum multiplied by the transport 
cross section, enters into the direct terms. and I n , which involve the t-matrix elements 
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themselves and not just their absolute square, appear only in the exchange terms, i.e., 
those proportional to e, which arise from (anti)symmetrization of the wave functions. To 
study the relative importance of the latter it is necessary to specify the 2-body potential 
energy function for the problem of interest and calculate the i-matrix that it produces. 
For He the interatomic potential becomes strongly repulsive at distances less than 2.4 A 
[fiH l, so we shall use a hard sphere potential as a guide; and there are two limiting cases 
for which analytic approximations to the scattering amplitude can be obtained: ka/fi 1 
and ka/h <C 1 where k is the relative momentum of the colliding particles and a is the 
radius of the hard sphere potential. Numerical calculations of the quantum mechanical 
scattering amplitude in the short wavelength limit [15] are in qualitative and even semi- 
quantitative agreement with the classical scattering of a scalar wave on a rigid sphere, as 
presented in Morse and Feshbach [[16||, and also by the scattering of an electromagnetic 
wave by a perfectly conducting sphere, given in Jackson [17]. That t-matrix takes the 
approximate form 

. 2vrn 2 a fka (1 + cos 0) Ji(¥ sin0) ka 



t 



a / na (i -t- cosciJii-*- amu . Ka a . \ _ , .„„, 

"oVlT hT~~lT - + exp-2i— sin- ) (ka h » 1 37 



where the first term arises from scattering by the shadow side and the second term from the 
illuminated side; /x is the reduced mass. The pure imaginary shadow scattering amplitude 

dominates in the small angle region 9 ~ h/ka, and oscillates rapidly with decreasing 
amplitude for larger angles. 



Normalizing once again to the conditions in reference [12] gives 

^ = 12 . 5 (J!LJ1)" 2 « (38 , 
h VM3 He 293 / 2.4A V ; 

where we have chosen 2.4 A as the hard sphere radius for a gas of 3 He atoms. Although 
it would require a very low temperature to get into the long wavelength limit, even a 
moderately low temperature, ~ 10°, would introduce significant corrections to the i-matrix 
as well as the first order collision terms, S v ^\ In particular, contributions to S v ^ arising 
from antisymmetrization of the wave functions would not be negligible and this would 
produce important modifications to the remainder of the paper. Keeping that in mind 
we now estimate the magnitudes of Iu,Ii and I n in Eq.(A-23) using the short wavelength 
expression (|37|). It should be noted that this limit, A/a <C 1, together with the low density 
assumption, na 3 <C 1, imply nX <C 1, which is the negligible degeneracy assumption made 
earlier. Using Eq.(|37[) we can also verify the assumption made earlier that the shift in 
the momentum (or momentum transfer) argument of the t-matrix element by the amount 
jhF(t) has a negligible effect for the conditions of reference |jp2]| . 

Although the shadow scattering and the scattering from the illuminated side each 
contribute ira 2 to the total cross section, the factor (1 —cos 9) in Eq.( A-22]) for the transport 
cross section eliminates the shadow contribution, hence o~u(k) = na 2 ; and from Eq.( [A-23| ) 
Ijj = Tik^a 2 , where A: 3 is the average of A; 3 over a Maxwell-Boltzmann distribution for a 
particle having mass [i. It takes a considerable amount of algebra to show that when the 
short wavelength expression Eq.(37) for the scattering amplitude is inserted into Eqs.([A^| 
Jj), (|A^22|) and flA^f ), // is smaller than I v by 0[(ka/h)- 2 }. I n is even smaller, being 
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0(exp[— (ka/h)~ 2 ]). Therefore all the terms proportional to e in Eqs.( A724| )-( |A-26| ) can 
be neglected. 

In the short wavelength limit Eqs. ( A-24 )-( A-26 ) for the time derivative of the first 
order functions h^ 1 ' become 

n 1/2 ^ = n 1/2 h(3^G(t) + ^^n 1/2 n_ 1/2 [h% 2 - hf/ 2 ]I v (39) 
n-i/2—^ 11 = -n_ 1/2 h(3^G(t) - ^^n 1/2 n_ 1/2 [h W 1/2 - h$ 2 ]I u (40) 



and 



= l Jd m + W^ 2 ^ + n ~W h % ~ ("1/2 + n- 1/2 )h^\l v . (41) 

Summing the first two equations gives 

n l/2 h { ^ 2 + n. 1/2 h^\ /2 = (n 1/2 - n_ l/2 )h(3^ F{t) (42) 

where F(t) is the time integral of the field gradient, Eq.([i~9|), There is no need to solve 
4/2 an d ^-1/2 



for /i^ 1 ! and hPJt /2 separately because only the combination shown in Eq.(|42|) is needed in 



the equation for dhft' /dt; and furthermore the terms of first order in 7 do not contribute 
to the magnetization. 

The solutions for the real and imaginary parts of h+ are 



Reh + } = ^W" (ni/2 " n -i/2)^ (M)/ c/ exp(-ot) J* dt'F{t') exp(at') (43) 

and 



1 r* 

Im/jW = _ exp(-at) J dt'F(t')exp(at') (44) 

where 

4/3 (M) 

a = — Ti—nlu- (45) 
3M v ; 

n = n\/ 2 + n_ij 2 is the total number of particles per unit volume. 

a can also be expressed in terms of the diffusion constant, D, for a related problem. 
Suppose there is no magnetic field and no transverse component of the magnetization, 
but there is a spatial gradient to the z-component of magnetization. This reduces to a 
standard 2-component problem with equal masses (spin up and spin down); and if the 
diffusion constant is calculated with the same approximation as above, namely keeping 
only the leading term in an expansion in Sonine polynomials, then [18] 
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The qualitative significance of or 1 can be seen from Eq.(f45|) or its expression in terms 
of the diffusion constant: it is the mean free time between collisions. Normalising to the 
conditions in reference IFC 



a 



Provided G(t), and therefore F(t), vary on a time scale much longer than a -1 , the integrals 
in Eqs.(|43|) and ( |44| ) are well approximated for t » a -1 by 

exp(-at) [ dt'F(t')exp(at') = — . (48) 
jo ol 

This leads to 

Re/i (i) = n 1/2 -n. 1/2 ^ (M)F 
n 

and 

Im/jW = ^f(t) (50) 
+ Ma w y ' 



For conditions similar to those of reference |l2l] Re/i^ is much smaller than Im/i^, their 



ratio being, at most, that of the mean wavelength to the mean free path; this will be small 
if the previous assumptions that Ac« and a <C d are both satisfied, and we shall neglect 
Re/i^ . It vanishes exactly if there is no z-component to the expectation value of the spin. 

Second Order 

There is a restriction on the form of a distortion of the distribution functions f±\/2 that is 
of second order in 7 (and therefore also of second order in the direction vector u), which 
arises from the requirement that it not alter the values of the particle number density, n, 
and the expectation value of the z-component of spin, < a z >, both of which are fixed by 
the terms of zero order. Since these quantities are obtained from Eqs.(|i~5|) and (|i~6|), the 
integral over all momenta of the distortions must vanish. The only allowed form is 

/g /2 (p) = n ±1/2 B( M Hp)[u • u - 20M(u • P) 2 ]/4V (51) 
Indeed, such a form arises naturally from a second application of the operator u • V p from 
!4) on f±y 2 - ^ or 9+^ we wr ite 

gf =n + B( M \p)(u-p) 2 h^(t), (52) 



which arises similarly from a second application of u • p from Eq.(^6[) on g+ 



Now pick out of Eqs.(|24|)-(|26|) for the time derivatives of the distribution functions the 

(2) 

terms of second order in 7. When the equations for df±y 2 /dt are integrated over all p 

(2) 

they become 0=0. The equation for the time derivative of g + becomes 

,,(2) 

n +J BW(p)(u • P ) 2 ^ = i±n + B( M )(p)(u ■ p) 2 F(t)h? + S V + (p) (53) 
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Integrating this equation over all p and changing the variables of integration in the collision 
term to (P,k) that term can be shown to contribute nothing; hence 



. i 

'M 



dtF{t)h { l\t). 



Inserting the expression for Ivsxh^i) from Eq.(50) gives 
1 rt 



M 2 a Jo 



dtF\t). 



(54) 



(55) 



4 Time Dependence of < cr > 

We now gather together all the contributions to < a >, the expectation value of the 
particle's spin, up to second order in the strength of the gradient of the magnetic field, by 
inserting into equations (p|) and ( 17| ) the results from fl2C|), fl27| ) -(|3l|) , (p^) and (|55|). 



and 



< a z >=< a z > (t = 0) 



< 0+ >=< cr + > (t = 0) • exp[-i~/(B Q t + u • xF(i))] 



(56) 



1-7' 



l D f dtF 2 (t)+0( 7 4 ) 
Jo 



(57) 



where we have made use of the expression for the diffusion constant D from Eq,(|46|). 
Eq .(|57|) represents the rate of decrease of the expectation value of the transverse component 
of the particle's spin to second order in the strength of the gradient of the magnetic field. 

How does this solution for < a > compare with the procedure followed in references 
[|J-[|J] mentioned at the beginning of the paper? This consists of adding a classical diffusion 
term to the equation for the time rate of change of the magnetization, i.e., 



dM 
~di~ 



7(M x B) + DV 2 M. 



(58) 



Recall that M is given by < a > multiplied by the magnetic moment of a particle, "fh/2, 
and the number of particles per unit volume. With B given in Eq.(^) and the initial value 
of M independent of postion, the exact solution of Eq.([5^) is that M z remains at its initial 
value; and M + is given by 



M+(x, t) = M + (t = 0) exp[-i~f(B t + u • xF(t))] exp 



7 2 L> /* F 2 {t)dt 
Jo 



(59) 



The final exponential factor is found in references 0] and [12]. If that exponential is 
expanded for small argument, to second order in 7 it gives the corresponding term from 

It was stated in the Introduction that the rate of decrease of the transverse component 
of magnetization is affected by both thermal motion and the interaction of a particle 
with its environment. By turning off the 2-body potential in Eq,(p6|), thereby setting the 
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collision terms S to zero, the rate due just to the random velocity distribution can be 
found, and the exact result is 



M + (x, t) = M+(t = 0) exp[-i~f(B t + u • xF(t))] exp 



,kT 



t \ 2 n 

F(t)dt 



M Wo 

(no collisions) (60) 



Using expression ( jig ) for the diffusion constant shows that the ratio of the arguments of 
the exponentials in Eqs.fl6C|) and (]59| ) is qualitatively given by aT where a -1 is the mean 
free time between collisions and T is the duration of the gradient. For the conditions 
of the experiment in reference [12] aT ~ 10 7 . It is not surprising that the transverse 
component of magnetization decreases more rapidly without collisions because the particles 
sample a much larger range of magnetic field strengths in a given time than they would if 
experiencing collisions, and consequently their precession rates cover a larger range. 

For the conditions of the experiment in reference [12] the quantity Jq F 2 (t)dt has 
the value 0.032 sec/mm 2 , where T is the time at which the gradient is turned off. Using the 



experimental values 15] of the ratio |M + (T)/M + (0)| of the magnitudes of the transverse 
magnetization (at various points along the length of the cylinder) after the gradient is 
turned off to its initial value, and comparing with Eq. (|59|) , the value of the 3 He self 



diffusion coefficient is found to be D=15.9 mm /sec [15|. The discrepancy with the values 



shown in Fig. 5a of reference [12] is due to a numerical error in that reference [15|. Although 
the statistical error in the value of D is small, only 0.4 mm 2 /sec, there is a much larger 
systematic error of almost 25 percent because the cell was prepared initially at much higher 
than room temperature |l5| ]. 

We also evaluated D from its relation to a, Eq. (f46|) , using the estimate given in (f47|) 
based on a hard sphere potential with radius a=2.4 A, and find D= 23.1 mm 2 /sec. This 
falls between the results obtained in previous experiments [fiT 



and [20 1, after correcting 
for pressure and isotopic mass dependence, and that reported in |l2|, O]. 



5 Summary 

The system under consideration consists of spin- 1/2 particles in a magnetic field with 
a time-varying gradient to its strength, and with its initial magnetization, taken to be 
spatially uniform, having components both transverse to and along the field direction. 
We now summarize the line of reasoning that led to the expression (^?j) for the time 
development of the expectation value of the transverse component of spin, < cr+ >, starting 
with Eq.(||) for the time derivative of p, the 1-body density matrix. It is the off-diagonal 
elements of p in the 2-dimensional spin space, or equivalently the Wigner distribution 
function / + , that determine < a + >. 

A number of parameters play a role in this development and certain inequalities among 
them were assumed, in succession. As each of these relations arose we compared the rele- 
vant ratios with their values in the experiment reported by Schmidt et al. fll2|] ) expressed 
in a way that allows quick comparison with any other set of experimental parameters. For 
a low density system of spin 1/2 particles, in the equation for the 2-body density matrix 
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we neglected 3-body collisions, and this requires that the mean spacing between particles, 
d, be large compare to the range of the force, a. We also assumed that prior to a collision 
the particles are uncorrelated, so that the 2-body density matrix is the product of the 
respective 1-body density matrices. Quantum scattering theory is then used to find the 
effect of a collision. 

After expressing the 1-body density matrix elements in terms of Wigner distribution 
functions, / , which form a 2x2 matrix in spin space, a set of coupled nonlinear integro- 
partial differential equations ([ji]) for the /'s is obtained. On the right hand side are three 
contributions to their time rate of change: S arising from the kinetic energy; S" 7 from the 
magnetic energy; and S v from the 2-body potential energy terms in the Hamiltonian. The 
first two terms are simple to write down and are given in Eqs.(p|) -(|i"4]); the expressions 
for the collision terms S v are obtained in the Appendix, Eqs.( |A-17[ ) and QA-18 ). Two 
approximations were made to obtain those collision terms. The first is to assume that the 
2-body potential is purely central, and this is satisfied extremely well for 3 He atoms. The 
second assumes that the gradient of the magnetic field strength is sufficiently small; this 
arises from the spatial dependence of the phase factor that governs the Larmor precession, 
Eq.(20), where F(t) is the time integral of field gradient. We expected and verified that 
this is the only spatial dependence of the distribution functions. The product jF(t) plays 
the role of a wave number, and leads to the off diagonal element (in spin space) of the 
density matrix being non-local in momentum space by the amount jhF(t); see Eq.(p^). 
This nonlocality shifts the arguments of g+ and the i-matrix elements in the collision terms 
by amounts of the same order, and even leads to non-energy-conserving matrix elements. 
If these shifts were sizable they would represent important quantum mechanical effects. 
The magnitude of the effect of the shift on g + depends on the ratio of jhF(t) to the mean 
thermal momenta, and in Eq.(23) this is shown to be very small for the conditions of the 
experiment in reference [12]. Using a hard sphere potential as a guide, this same small 
ratio makes the effect of the shift in the arguments of the i-matrix elements negligible as 
well. 

To solve the resulting equations (p4|) - (26) we expanded the distribution functions 
in a perturbation series to second order in the strength of the field gradient, G(t), using 
the gyromagnetic ratio 7 as the expansion coefficient since it multiplies G; this makes 
the equations linear in each order. The only distortions away from Maxwell-Boltzmann 
distributions that were used are of the form u • p in first order and u • u and (u • p) 2 in 
second order, where u is the direction of the gradient. Omitting additional dependence on 
p 2 is equivalent to keeping only the leading term in an expansion in Sonine polynomials 

& 

To first order in 7 the equations ( |A-24|) - ( |A-26|) were written still containing the effect 
of antisymmetization of the wave functions on the collision terms; these are the terms 
containing e and they involve, through the quantities and I n defined in ( [A-23 ) , not 
just the absolute square of the t-matrix elements but also their real and imaginary parts. 
For a hard sphere potential in the short wavelength limit, however, the exchange terms 
are smaller than the direct terms by 0(X/a) 2 , where a is the hard sphere radius, and we 
expect this to be true for a more general potential with a the range of the force. But note 
from Eq.(^) that for 3 He at room temperature the ratio a/A is only 12.5 with a = 2.4 A. 
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At a significantly lower temperature corrections due to the exchange terms would become 
important; they were neglected in the remainder of the paper. The direct terms, on the 
other hand, manifest themselves through the single quantity Ijj, which is the mean of the 
product of the cube of the momentum with the transport cross section ajj, averaged over 
a Maxwell-Boltzmann distribution for a particle having the reduced mass. 

Provided the field gradient G(t) varies on a time scale much longer than the mean free 
time between collisions, a -1 , the solution for the first order distortion of the transverse 



Wigner distribution function is given in Eqs.(49) and (|50|). ct is simply related to the self 
diffusion constant D for 3 He, which can be more readily realized in a different experiment 
with no field gradient but rather a spatial gradient to the z-component of magnetization; 
this relation is given in Eq.(fl6|). The solution for the second order contribution to the 
transverse distribution function follows readily from the first order term, and is given in 
Eq.(H). 

If all the restrictions mentioned above are satisfied, then the expression for the rate 
of decrease of the transverse component of spin, to second order in the strength of the 
gradient, is given in Eq.(57). Comparison with the classical treatment described in the 
Introduction, the exact solution for which is given in Eq.(^), shows agreement to second 
order in the strength of the gradient. It was also shown that the rate would be much 
greater if there were no collisions, because the particles would sample a much larger range 
of magnetic field strengths in a given time and consequently their precession rates would 
also. 

It is not possible to make a good comparison between the magnitude of the 3 He self 
diffusion coefficient estimated from a hard sphere potential using Eqs.(H3) and ([47]), and 



the experimental result in [12] because the temperature varied over a large range [15|. 
The theoretical estimate of D= 23.1 mm 2 /sec falls between the values obtained in that 
experiment and that obtained in previous experiments [19, 20 j. 
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A Appendix 

Here we work out the expressions for the collision terms S v . After substituting Eq. 
p( 2 )(l,2) in Eq.(|l~8l) the elastic scattering t-matrix is introduced via the relations Jll 
fi( + ) = (I + G^ + H) and Vfl^ = i, where is the Green's function. Since the potential 
conserves the total momentum and the i-matrix elements depend only on the relative 
momenta of the two particles, they can be written 

< p A ,m A ;p B ,m B \t\p' A ,m' A ;p' B ,m' B > 

■(2irhf5^(p A + p B - (p' A + p' B )) t(*AB, k'AB) (A-l) 
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where ]s.ab = (pa ~ Pb)/2. The corresponding equation for the matrix element of is 

< p A ,mA;PB,rn B \0,^\p A ,m' A ;p' B ,m' B > = 5 

.(2nh) 3 5®(p A + PB- (p' A + p' B m {+) ^AB, k' AB ) (A-2) 

where 

rt+\k AB ,k AB ) = [(2vr/i) 3 <5 (3) (PA - p' A ) + t(k AB ,k' AB )/(E' AB +ir- E AB )\ (A-3) 

with Eab the sum of the 1-body energies of particles A and B. (It is understood that the 
limit r — > is to be taken.) Since the magnetic energy term in is the same in the 
primed and unprimed states only the kinetic energy contributes to the denominator in 
Eq. flA^I) . Hence E A b 

can be taken to be k AB /2fi where (j,, the reduced mass, is M/2. 
As described in the text we are assuming that the magnitude of r yhF(t) is sufficiently 
small that the shift in the argument of the 5- function in Eq.(^) can be neglected. Hence 
we shall write 

< p,m\p\p',m' >= (2TThfS^(p - p')/w(p) (A-4) 

where 

p mm (p) = (2vr/i) 3 / m (p) (A-5) 

and 

P-i/2 1/2 (P) = (27rhfg + (p) (A-6) 
When two sets of intermediate states are introduced into Eq. (|l8|) two of the integrations 



can be carried out immediately by making use of Eqs. flA- H ), (A-2), and (A-4), giving 

1 sr^ f d 3 p2 d 3 p' d 3 p' 2 



5 ™ m ' (x ' p) ~ h(27rh) 3 ^I 



h (2nh) 3 f^J (2nh) 3 {2irh) 3 (2nh) 
•(2^)V 3 )(p + p 2 -(p' + p 2 )) 

{[t(k,k')^ (+) *(k,k') -n(+)(k,k')t*(k,k')] Pmm>(p')Pm 2 mM) 

+e[t(-k,k')^ (+) *(k,k') - OW(-k,k')f(k,k')] p m2 m>{p')Pmm 2 (P2)} (A-7) 

After substituting for from Eq. ( |A-3| ) , for the terms linear in the f-matrix elements 
the integrations over the primed variables can be carried out, and for the terms bilinear 
in t the energy denominators can be combined, giving 

5mm' ( X >P) = Y\ (2TTh) 3 ^ J {2lTh) 3 V^' ^ ~ t *( k ^ k )]Pmm'(p)pm 2 m 2 {P2) 

+e[t(-k,k)p (p)Pmm 2 (P2) - t*(-k,k)p (P2)Pmm 2 (P 
1 y-v f d 3 p 2 d 3 p' d 3 p' 2 



{2nhf^J (2ttTi) 3 (2iTh) 3 (2-kTi) 3 



■(2TTh) 3 5^(p + P2 - (p' + p' 2 ))2m5(E' - E) 

• [|t(k, k') 1 2 p mm > (p>m 2 m 2 (p' 2 ) + (-k, k')t* (k, k')p OT2 m' (p')pmm 2 (p' 2 )] j (A-8) 
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This establishes the fact that only energy conserving matrix elements contribute to the 
collision terms, so i(k, k') can be written as t(k,9). 
Note that any integral of the form 



f d 3 p f d 3 p' 2 

■(27rft) 3 ^ 3 )(p + P2 - (p' + p 2 ))27n5(£' - E)F(p, p 2 , p', p' 2 ) 



(A-9) 



can be rewritten by changing the variables of integration to (P',k'), where P' = p' + p 2 
and k' = (p' — p' 2 )/2. After carrying out the integrations over P' and the magnitude of 
k', using E' — E = (k' 2 - k 2 )/2/j,, / becomes 



/(p ' P2) = 27r /{2W F(P ' P2 '^ 



(A-10) 



where Q,' stands for two angles specifying the direction of k'. The combination fikdQ' / {2irfi) i 
is the density of final states. 

Making use of the relation between the differential scattering cross section and the 
t-matrix elements gives 



= — |t(k,k 



dou 


2ir 


dSV 




da A 






~ T 



r\\2 



fj,k 



(2vr/i) ; 



and 



dai 2ir 



[t(k,k')+et(k,-k')]/V2 



fxk 



(2nhY 



T^- k '^ k ' k 'W 



(A-ll) 
(A-12) 

(A-13) 



where ajj and a a are, respectively, the unsymmetrized and antisymmetrized cross sections. 
a i is not actually a cross section, just a convenient notation for the interference term; 
v = k/fi is the relative velocity of the colliding particles. 

Further progress is made by making use of rotational invariance of the t-matrix, and 
unitarity of the S'-matrix, which takes the form 

■(2wh) 3 S< 3 \p + p 2 - (p' + p' 2 ))2mS(£' - E), (A-14) 
and from the expressions above this becomes 

2i Imi(k,k) = -iti I dn'v^j-. (A-15) 

J til L 

The corresponding equation for t(k, — k) is 

2i Imt(k, -k) = -iti I dQ'v-^j. (A-16) 

J CO, L 



18 



In terms of these quantities the separate cases of S v from Eq.( A-8| ) can be obtained by 
carrying out the summation over m 2 and rearranging some terms. 

SYptp) = J d 3 P2 dQ'{[f l/2 (p')f_ l/2 (p' 2 ) - /i /2 (p)/-i/ 2 (p 2 )]^ 

+[/l/2(p')/l/2(P 2 ) - /i/2(p)/i/ 2 (P2)]^^7 + es + (p')0+(p 2 >^7} 

-e^r J d 3 p 2 [t(k, -k)g + (p)g* + (p 2 ) - c.c] (A-17) 

The corresponding equation for S^ 1 ^ 2 ^ s obtained from ( |A-17 ) by making the interchanges 
fi/2 ^ /-1/2 an d g+ «-> g* + - Finally, the expression for becomes 

Sl(p) = expMBot + u • xF(t))] S£(x, p) = / rfW^' 

•{(5 + (p')[/l/2( P2 ) + /-l/ 2 (p 2 )] - 5 + (p)[/l/2(P2) + /-1/ 2 (P2)])- </ "' 



+eb+(p')/-i/2(p 2 ) + /i/2(p')9+(P2)]w ^7} 

i 



| d 3 p 2 (t(k,-k)[ 5+ (p)/_ 1/2 (p 2 ) +/i /2 (p) 5 +(P2)] 
t*(k, -k)[ 5+ (p 2 )/_ 1/2 (p) + /i/ 2 ( P2 ) 5 +(p)]) (A-18) 



We now extract the collision terms of first order in 7, to be used in Eqs.(|34|) - (p6|), by 
inserting Eqs.(|32|) and (33) in the above. It is helpful to express each individual momentum 
variable in terms of the total and relative momenta P,k and k'. When this is done the 
terms in S^w^ i nv °l vm g daA/dQ,' vanish by conservation of momentum; and making use 

of unitarity, Eq.( A-16 ), the terms involving the real part of h^P cancel, leaving 

s^\ P ) = U .J rfV!J'#>(p)BWwK«-i /2 (p' - p)(^ 1 / ) 2 " h% 2 )v^r 

+e\n + \ 2 2ik'lmh^v^} 

+e|n + | 2 ^Im/i ( h 1) u- J d 3 p 2 B^ M \p)B (M \p 2 )kRet(k, -k). (A-19) 
Furthermore, »S^ 2 (p) = ~~ ^1/2 (p)- Following similar reasoning leads to 



5^ (1) (p) = u - n + J d 3 p 2 dn'B ( - M \p)B ( - M \p 2 ) 

•{(p' - P)[^+ ) (wi/ 2 + n_i /2 ) - ("1/2^1/2 +n_i/ 2 '»?i/2)]«^' 
+ek'[(n 1/2 h { ^ 2 - n_ l/2 h {l) l/2 ) - (n 1/2 - n -i/2) /l + ) ] w ^} 
-ie-n + [{n 1/2 h^J 2 - n_ 1/2 h W 1/2 ) - (n 1/2 - n_ 1/2 )h { ^} 

■u- J d 3 p 2 B (M \p)B^(p 2 )kRet{k,-k) (A-20) 
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When these expressions for the first order collision terms are inserted into Eqs.(£ 
fl36|) we see that every term is of the form u- some vector, where u is an arbitrary direction 
in space; therefore those equations must be valid as vector equations with the u omitted. 
A necessary consequence can be obtained by forming p- those equations and integrating 
over all p. At the same time we make the following observations about the collision 
terms. Change the variables of integration from (p, P2) to (P, k), where P = p + P2 and 
k = (p — P2)/2, and make use of 

B {M) (p)£(M) = B {2M) (p)£(M) ( fc ) ( A _2i) 

which follows from conservation of energy. For each fixed value of k choose its direction as 
the polar axis for the integration over the angles of k', so that dO,' becomes d£l = 2ir sin ( 
where 9 is the angle between k and k'. We also define the following quantities. 



au {k) = [ cft2(l-cos6>) 
Jo 



day 
dfl 



ai(k) = j dfi cos 8^- (A-22) 



and 



I v = I d 3 kk 3 B^(k)au(k) =< k 3 a v (k) >W 
Ij = J d 3 kk 3 B^{k)a I (k) =< k 3 aj(k) >M 

J w = / d 3 kk 2 B (ll) {k)Ret(k,9 = vr) (A-23) 



From the definition of daj/dQ, Eq.( |A-13 ), it is seen that cri(k) and ij are pure imaginary, 



whereas Ijj and /„- are real. au{k) is called the transport cross section [ 13 J and Ijj is the 
mean value of k 3 au{k) averaged over a Maxwell-Boltzmann distribution for a particle of 
mass [A. 

Eqs.(|34|) - ([?(]) now become 

ni/2-^ = n l/2 n(i^G(t) 



+ 



dt 
3M 



[n l / 2 n_ l / 2 [h ( il /2 - hf/ 2 ]Iu - 2e|n + | 2 Im/i^_ 1) Im//| 



+£ 3^~ |n+l m/l + (A " 24) 
= -n_ 1/2 hf3^G(t) 

4/3( M ) 
3M 

-e^^\n + \ 2 lmh^ } I n (A-25) 

Oil 



[n 1/2 n_ 1/2 [h^l /2 - hfj 2 ]Iu - 2e|n+| 2 Im/i^ ) Im//| 
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and 

,(!) , „ . ..1,(1) . 



dhV 1 4/?( M ) r 

~dT = l M F ^ + IF U ni /2^V2 + "-1/2^1/2 - ("1/2 + n-l/2)>l+ ; Ktf 



+e[(ni/2^i/ 2 - ^-1/2^^1/2) - ( n i/2 - "--i/2)^+V/} 

-it 3h [("1/2^1/2 - «-l/2^_i/ 2 ) - ( n l/2 - "-1/2^+ Rtf- (A-26) 

These are four coupled ordinary linear differential equations for the four unknown functions 
^±1/2 an( i the real and imaginary parts of h+\ 
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